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Abstract: 

This  year’s  work  was  concerned  with  the  electronic  structure  of  ZnO  and  related  alloys. 
The  nature  of  the  valence  band  fine  structure  was  studied  and  progress  made  in  identifying 
the  problems  underlying  its  identification  by  experimental  optical  studies.  Valence  band 
effective  mass  parameters  were  determined.  The  electronic  structure  of  MgjjZni-j,©  alloys 
were  studied,  their  bandgap  bowing  determined  and  their  band  offsets  to  ZnO  and  other 
related  materials  were  investigated.  The  stability  of  MgO  in  the  wurtzite  structure  relative 
to  the  rocksalt  structure  was  investigated.  A  continuous  transformation  from  rocksalt  to 
wurtzite  was  discovered  and  applied  both  to  MgO  and  to  GaN. 
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1  Introduction 


This  report  describes  research  performed  under  grant  no.  NOOO 14-98- 1-01 60  during  the  performace 
period  Nov.  1998-Sept.  1999. 

Within  the  general  topic  of  modeling  of  alloys  of  wide-band-gap  semiconductors,  this  year’s  work 
was  focused  on  the  wide  band  gap  semiconductor  ZnO  and  its  related  alloys  with  Mg,  which  have 
recently  attracted  considerable  attention.  The  reasons  for  this  renewed  interest  in  ZnO  are:  1)  the 
increasing  availability  of  large  and  high-quality  bulk  single  crystals  of  ZnO  by  a  variety  of  growth  meth¬ 
ods,  including  low-cost  hydrothermal  growth;  2)  the  recent  demonstration  of  optically  pumped  stim¬ 
ulated  emission,  [1]  3)  the  development  of  epitaxial  growth  methods  of  ZnO  and  related  Mga;Zni_a;0 
alloys,  [2]  4)  the  realization  that  ZnO  is  closely  lattice  mactched  to  GaN  but  may  be  an  interesting 
electrooptical  material  in  its  own  right,  surpassing  GaN  in  certain  aspects.  The  first  International 
Workshop  on  ZnO  was  held  in  Dayton,  Ohio,  from  October  7-8,  1999.  The  high  quality  of  the  crystals 
of  ZnO  is  demonstrated  among  others  by  the  possibility  to  obtain  fine  structure  information  on  the 
valence  band  structure  using  photoluminescence.  [3]  This  work  has  re-opened  a  30-year  old  controversy 
over  the  ordering  of  the  states  near  the  valence  band  maximum.  Our  work  described  below  was  in 
part  geared  towards  a  better  understanding  of  this  controversy. 

Our  work  has  focused  on  two  topics:  1)  the  valence  band  ordering  and  2)  the  properties  and 
stability  of  the  alloys  with  MgO.  Important  progress  has  been  made  on  both  aspects  and  been  presented 
already  at  some  meetings.  Some  further  work  is  necessary  to  bring  these  investigations  to  a  definitive 
conclusion.  Thus,  a  nocost  extension  was  requested  to  continue  this  work. 

We  now  briefly  describe  the  progress  made  under  each  of  these  topics  and  then  provide  lists  of 
presentations  and  publications  obtained  or  in  progress  under  this  work. 


2  Valence  band  structure 

We  here  describe  the  problem  of  the  valence  band  structure  along  with  our  specific  calculations  which 
provide  progress  towards  its  resolution. 

The  nature  of  the  valence  band  maximum  fine  structure  in  ZnO  was  first  studied  by  Thomas.  [4] 
Based  on  the  polarization  dependence  of  the  absorption  and  reflectivity  spectra,  he  came  to  the 
conclusion  that  the  valence  band  ordering  of  ZnO  is  inverted  from  the  usual  one  in  other  II- VI  wurtzite 
materials.  Namely,  the  symmetry  character  of  the  highest  valence  band  state  (leading  to  the  so-called 
A-exciton)  according  to  his  analysis  was  Tt  rather  than  Tg.  The  symmetry  of  the  B  state  would  then 
be  Tg.  The  C-exciton  also  has  symmetry  Tt.  These  splittings  come  about  as  a  result  of  the  interplay 
between  spin-orbit  coupling  and  crystal  field  splitting.  His  interpretation  was  that  the  crystal  field 
splitting  is  large  compared  to  the  spin-orbit  splitting.  Without  spin-orbit  coupling,  the  valence  band 
would  be  split  into  a  singlet  Ti  and  a  doublet  Ts.  The  singlet  Ti  is  essentially  a  pz  like  state  (with 
slight  s-admixture)  while  the  Ts  is  a  pi,  Py  like  state.  With  the  usual  Fs  >  Fi  ordering,  one  thus 
expects  the  C-exciton  to  be  allowed  in  E  ||  c  polarization  and  the  F5  derived  (unresolved  A,B  excitons) 
to  be  allowed  for  the  other  polarization  E  J.  c.  This  is  essentially  what  was  observed  in  Thomas’s 
reflectivity  spectra.  Absorption  reveals  the  fine  structure  splitting  in  A  and  B  excitons.  With  the 
spin-orbit  coupling  present,  the  F7  state  derived  from  the  F5  will  oi'tain  a  slight  admixture  of  while 
the  Fg  stays  purely  x  and  y  like.  Thus,  one  expects  the  Fy  state?  to  b€eo;i3e  wsakly  oil  wed  for  E  ||  c. 
Thomas’s  data  revealed  this  to  be  the  case  for  the  highest  energj  A-excJton.  Bence,  liis  assignment. 
This  was  later  confirmed  also  by  by  Liang  and  Yoffe.[5] 
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On  the  other  hand,  the  interpretation  of  the  lines  by  Thomas  was  challenged  by  Park  et  al.  [6] 
They  concluded  on  the  basis  of  temperature  dependence  studies  and  on  the  basis  of  a  prior  study  by 
Reynolds  et  al.[7]  that  the  lines  which  were  identified  as  the  A-exciton  lines  Am  and  Al  (the  latter 
being  a  longitudinal  exciton)  by  Thomas,  were,  in  fact,  not  arising  from  a  fi:ee  exciton  but  corresponded 
to  the  Ps  and  Pe  states  of  a  ionized  donor  bound  exciton.  This  was  later  disputed  by  other  authors, 
among  them  B.  Segall,[8]  who  confirmed  the  free  exciton  nature  of  the  absorption  line  in  question  by 
a  study  of  the  phonon  assisted  absorption  onset.  The  controversy  was  never  fully  resolved,  although 
most  future  papers  adopted  the  interpretation  of  Thomas. 

Recently,  Reynolds  et  al.  [3]  were  able  to  restudy  this  issue  using  second  order  photoluminescence 
spectra.  The  second  order  of  the  diffraction  grating  provided  them  higher  resolution,  fully  resolving 
the  additional  fine  structure  of  the  excitons  due  to  the  exchange  splitting.  Furthermore,  they  studied 
the  behavior  of  these  spectra  in  a  magnetic  field.  Their  conclusion  was  that  the  line  is  indeed  a  firee 
exciton  line,  but  that  nevertheless,  the  valence  band  maximum  must  have  Pg  symmetry.  We  will 
discuss  their  arguments  below  in  more  detail. 

The  meaning  of  the  ordering  proposed  by  Thomas  is  that  the  spin-orbit  splitting  is  negative.  The 
possibility  of  a  negative  spin-orbit  splitting  was  first  suggested  by  Cardona  in  a  study  of  copper  and 
silver  halides.  [9]  It  was  subsequently  explained  by  Shindo.[10]  The  origin  is  the  presence  of  lower  lying 
d-bands.  The  valence  band  maximum  being  an  antibonding  combination  of  anion  p-like  states  and 
cation  d-like  states,  results  in  a  negative  contribution  of  the  atomic  d-orbitals  to  the  effective  spin-orbit 
splitting.  Thus,  one  expects  the  possibility  of  a  negative  spin-orbit  parameter  if  the  d-bands  lie  fairly 
close  to  the  valence  band  maximum  and  have  a  strong  atomic  spin-orbit  parameter.  This  is  clearly 
the  case  in  Cu  compounds.  The  situation  is  marginal  in  ZnO,  because  the  d-bands  here  lie  about 
7  eV  below  the  valence  band  maximum  (according  to  photoemission  data[ll]).  In  previous  work  on 
nitrides,  we  already  pointed  out  the  importance  of  the  d-bands  in  reducing  the  values  of  the  spin-orbit 
splittings.  [12] 

Prom  the  above,  it  is  clear  that  the  results  of  a  negative  spin-orbit  splitting  may  depend  crucially 
on  the  position  of  the  d-bands.  Thus,  we  paid  special  attention  to  this  feature  in  the  band  structures. 
Our  initial  band  structure  calculations,  using  the  local  density  approximation  (LDA)  obtain  a  d-band 
position  of  about  5  eV  below  the  VBM  in  contrast  to  the  experimental  value  of  6.95  eV.[ll]  The  sign 
of  the  spin-orbit  splitting  in  this  calculation  was  determined  in  two  ways.  First,  in  zincblende  ZnO, 
the  valence  band  splits  into  a  4-fold  state  of  symmetry  Fg  and  a  2-fold  state  of  symmetry  Ft.  Thus 
simple  inspection  of  the  degeneracy  of  the  eigenvalues  revealed  that  F7  was  lying  above  Fg,  meaning  a 
negative  spin-orbit  parameter.  Secondly,  inspection  of  the  wave  functions  in  the  wurtzite  case  revealed 
that  the  highest  valence  band  contained  a  and  s  components,  indicating  F7  symmetry  while  the 
second  state  had  pure  Px,  Py  (and  some  d-admixture)  but  absolutely  zero  pz  or  s  components.  So, 
this  means  that  in  LDA,  a  negative  spin-orbit  splitting,  confirming  Thomas’s  conclusion  is  predicted. 

However,  one  may  argue  that  the  LDA  calculation  places  the  Zn  3d  band  too  high  in  energy  and 
thus  overestimates  the  negative  component  of  the  spin-orbit  splitting.  We  thus  further  investigated 
the  behavior  of  the  splitting  as  a  function  of  d-band  position.  These  results  are  shown  in  Fig.  1 
and  Fig.  2.  What  we  find  is  that  both  the  crystal  field  and  the  spin-orbit  splitting  depend  strongly 
on  the  d-band  position.  In  fact,  both  decrease  monotonically  and  nearly  linearly  as  a  function  of 
increasing  d-band  binding  energy.  We  find  good  agreement  with  the  experimentally  deduced  Ea  —  Eb 
and  Eb  -  Ec  splittings  for  a  d-band  position  of  -6.25  eV.  At  the  experimental  d-band  position,  the 
crystal  field  splitting  is  strongly  underestimated.  How  to  explain  this? 

The  origins  of  discrepancy  between  LDA  and  photoemission  on  the  d-band  position  are  twofold. 
First,  there  is  the  fact  that  LDA  treats  exchange  and  correlation  in  an  orbital  independent  manner 
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Figure  1:  Valence  band  splittings  versus  d-band  position  in  wurtzite  ZnO 


and  Ag^at  different  d-level  for  ZnO 


Figure  2:  Crystal  field  (Ai  =  Acr  and  spin-orbit  splittings  as  function  of  Zn  3d  band  position:  fits 
derived  either  with  the  quasicubic  model  using  cubic  Ago  =  3A2  =  SAa  or  with  the  anisotropic  spin 
orbit  parameters  A2  and  A3  are  shown. 


4 


and  leads  to  a  self-interaction  error  because  of  imperfect  cancellation  of  the  Coulomb  and  exchange 
integrals,  in  particular  for  localized  states.  Secondly,  however,  in  the  photoemission  experiment,  a 
final  state  relaxation  effect  is  involved.  In  some  sense,  in  order  to  obtain  the  correct  valence  band 
splitting  we  need  to  include  the  first  but  not  the  second  effect. 

Another  way  of  phrasing  this  is  within  the  language  of  quasiparticle  theory.  The  energy  of  a  state  as 
determined  by  photoemission  is  the  quasiparticle  energy,  which  follows  from  an  equation  which  contains 
a  non-local  and  energy  dependent  self-energy  operator  instead  of  a  exchange  and  correlation  potential. 
This  self  energy  operator  in  the  GW  approximation  is  essentially  like  a  dynamically  screened  Hartree 
Fock  exchange  term.  The  energy  dependence  of  this  term  implies  that  there  is  strictly  no  hope  of 
obtaining  both  the  valence  band  maimum  and  the  d-band  states  from  the  diagonalization  of  one  energy 
independent  Hamiltonian  matrix.  On  the  other  hand,  let’s  assume  we  make  a  static  approximation 
to  the  screening  of  the  exchange.  In  that  case,  we  will  already  correct  the  self  interaction  effect  that 
LDA  gets  wrong  and  obtain  a  much  better  orbital  dependent  treatment  of  exchange.  Essentially,  this 
will  lead  to  a  downward  shift  of  the  d-band.  This  is  what  our  calculation  is  simulating  in  a  simple 
manner.  We  simply  shift  the  d-band  atomic  orbital  energy  or  diagonal  matrix  element  in  our  LMTO 
representation  of  the  Hamiltonian.  Clearly,  following  the  above  argument  we  should  not  expect  that 
we  get  the  d-band  and  VBM  correct  with  the  same  shift.  Thus,  since  we  are  here  interested  in  the 
VBM,  we  adopt  a  slightly  empirical  point  of  view  and  use  this  shift  as  an  adoptable  parameter.  The 
result  that  we  obtain  a  d-band  position  intermediate  between  the  LDA  value  and  the  experimental 
value  for  the  correct  VBM  splitting  is  completely  in  agreement  with  expectations  of  the  above  outlined 
theory.  Of  course,  a  further  ab-initio  justification  of  the  d-band  shfit  required  would  be  preferable 
and  will  be  investigated  in  future  work.  We  have  in  fact  hopes  that  this  can  be  accomplished  within 
the  frame  work  of  recently  proposed  screened  exchange  Generalized  Kohn  Sham  methods. [13]  For 
the  moment,  however,  let  us  be  satisfied  that  indeed  both  the  Ea  —  Eb  and  Eb  —  Ec  splittings  can 
simultaneously  be  fitted  with  one  empirical  d-band  shift  parameter. 

Now,  we  can  with  this  shift  examine  again  the  symmetry  of  the  VBM  states.  We  find,  we  still 
have  Fy  symmetry  for  the  VBM.  Also,  if  we  use  the  same  shift  of  the  d-band  in  zincblende,  we 
still  find  the  doublet  above  the  quartet.  Finally,  Fig.  2  shows  that  assuming  a  negative  spin  orbit 
splitting  at  this  d-band  position  leads  to  the  expected  nearly  linear  dependence  vs.  the  d-band 
position,  whereas  assuming  a  positive  splitting,  would  lead  to  a  nonmonotonic  behavior.  Clearly,  we 
can  easily  identify  the  point  where  the  spin-orbit  splitting  passes  through  zero  as  a  function  of  d-band 
position.  Furhermore,  we  note  that  if  we  shifted  the  d-band  to  the  point  where  the  spin-orbit  splitting 
becomes  positive,  the  crystal  field  splitting  would  be  strongly  underestimated.  Further  corrections 
beyond  LDA  would  be  expected  to  perhaps  reduce  this  splitting  but  not  to  increase  it.  In  fact,  the 
gap  is  still  underestimated  and  including  a  gap  shift  will  reduce  the  interaction  of  the  Fi  with  the 
conduction  band  and  thus  reduce  the  crystal  field  splitting.  This  effect  is  expected  to  be  small.  We 
note  finally  that  with  this  d-band  shift,  the  gap  correction  required  for  ZnO  is  similar  to  that  in  GaN. 
This  is  expected  because  the  two  materials  have  nearly  the  same  experimental  gap,  as  well  as  close 
correspondence  in  lattice  constant,  density  of  electrons,  dielectric  constant,  all  factors  which  influence 
the  GW  corrections. 

So,  in  summary,  there  is  no  doubt  that  the  calculation  leads  to  a  negative  spin-orbit  splitting  even 
in  a  model  that  goes  beyond  LDA  by  iacuding  the  expected  shifts  of  the  d-band  and  their  influence 
on  the  crystal  field,  spin-orbit  splittir  s  ar.a’  gap  cc,  ri-t  ciioxi. 

It  then  remains  to  discuss  the  disi .  epancy  with  the  recent  conclusions  by  Reynolds  et  al.  [3]  To 
that  end  we  must  first  realize  that  their  conclusion  is  indirect  because  it  is  based  on  the  behavior  of 
excitons  rather  than  of  the  valence  bands  directly.  That  is,  one  might  have  that  the  exciton  binding 
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energy  of  the  A,  B  and  C,  derived  excitons  is  sufficiently  different  that  the  ordering  of  the  exitons 
bands  differs  from  that  of  the  valence  bands.  In  fact,  one  may  not  even  strictly  speak  of  B,  and  C 
derived  excitons,  because  the  exciton  binding  energy  is  of  order  60  meV  while  the  A-B  spUtting  is  only 
9.5  meV.  This  problem  was  already  mentioned  in  the  theory  paper  by  Hopfield  [14]  accompanying  the 
Thomas  paper.  This  aspect  deserves  further  investigation.  A  theory  for  dealing  with  these  aspects  is 
available  through  Baldereschi  and  Lipari’s  work[15]  but  implementation  of  it  for  wurtzite  remains  to 
be  accomplished.  We  hope  to  address  this  in  future  work. 

Another  explanation  of  the  discrepancy  however  is  possible.  While  there  is  of  course  no  question 
about  the  validity  of  the  data  of  Reynolds  et  al.,  their  assigment  of  the  features  to  the  symmetry 
should  be  carefully  examined.  When  one  considers  the  excitons  formed  by  a  Fy  valence  band  and  a 
Fy  conduction  band,  three  possible  symmetries  are  obtained:  Fi,  F2  and  F5.  When,  considering  a  Fg 
valence  band,  the  possible  symmetries  of  the  exciton  are:  F5  and  Fg.  In  short: 

F7®F7  =  Fi©F20F5 

F9®F7  =  FgSFg. 

Now,  F5  is  allowed  strictly  speaking  only  for  E  -L  c  and  Fi  is  allowed  for  E  ||  c.  The  other  lines  are 
forbidden  but  may  still  be  weakly  visible  because  of  not  exact  k-conservation,  or  because  of  not  perfect 
alignment  of  the  beams  vs.  the  crystal  axis.  The  argument  of  Reynolds  et  al.[3]  is  primarily  based  on 
the  behavior  of  PL  lines  in  a  magnetic  field.  What  one  finds  is  a  low  energy  line  at  3.375  eV  which 
exhibits  a  strong  splitting  in  a  magnetic  field  and  a  higher  energy  line  at  3.373  eV  with  negligigle 
splitting.  They  assign  the  strongly  split  line  to  Fg  because  one  expects  for  Fg  a  g-factor  which  is 
the  sum  of  the  electron  and  hole  g-factors.  This  then  leads  to  the  Fg  assignment  of  the  A-exciton. 
Their  g-factors  obtained  are  in  agreement  with  those  extracted  from  shallowe  neutral  donor  bound 
excitons  for  the  holes  and  with  the  known  electron  g-factor.  It  is  somewhat  difficult  to  explain  in 
their  interpretation  why  the  forbidden  Fg  line  is  almost  as  strong  as  Fg  for  the  E  ||  c  polarization  but 
not  for  the  E  ±  c  polarization,  although  one  might  in  luminescence  invoke  a  thermal  Boltzman  factor 
favoring  the  low  energy  state.  At  2K  the  Boltzman  factor  exp  {-AE/kT)  for  AE  «  2  meV  is  about 
a  factor  12. 

On  the  other  hand,  there  is  no  reason  why  the  g-factor  for  the  bound  exciton  should  equal  that 
of  the  free  exciton,  and  in  fact,  the  g-factor  for  the  holes  could  well  be  negative.  A  negative  g-factor 
was  reported  previously  for  ZnO  by  Blattner  et  al.  [16]  and  was  recently  also  reported  for  the  F7  hole 
in  GaN  by  Campo  et  al.  [17].  In  that  case,  an  alternative  interpreation  of  the  strongly  split  line  in 
magnetic  field  is  that  these  represent  the  Fi  and  Fg  lines.  The  Fg  is  forbidden  without  magnetic  field 
but  becomes  allowed  in  a  magnetic  field  because  the  latter  has  the  symmetry  Fg.  The  two  states  then 
couple  and  repel  each  other  resulting  in  splitting.  In  that  case,  the  Fg  derived  from  two  F7 ’s  should 
split  as  the  sum  of  the  g-factors  while  the  Fi  -  Fg  splitting  will  behave  as  the  difference.  So,  if  the 
F7  hole  g-factor  is  negative,  the  spectrum  can  alternatively  be  interpreted  as  being  derived  from  a 
F7  VBM.  In  fact,  this  interpretation  of  the  Zeeman  splitting  was  already  suggested  by  Thomas.  His 
Zeeman  study  was  done  on  absorption  rather  than  PL  lines  and  with  less  resolution  but  the  essentials 
stay  the  same.  One  might  ask  with  this  interpretation  why  the  Fi  is  so  weak  compaxed  to  Fg  for 
E  II  c?  However,  it  should  be  kept  in  mind  that  this  state  still  is  primarily  a;,y-like  with  only  a  very 
small  0  admixture  because  the  crystal  field  splitting  is  much  larger  than  the  spin-orbit  splitting.  So, 
even  with  a  Boltzman  factor  advantage,  it  may  be  quite  weak. 

In  any  case,  our  conclusion  from  this  whole  analysis  is  that  a  crucial  piece  of  missing  information 
is  an  indepenent  determination  of  the  sign  of  the  g-factor.  Thus,  we  have  initiated  work  to  calculate 
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Table  1:  Rashba-Sheika-Pikus  parameters*  for  ZnO  valence  band  effective  Hamiltonian 
Ai  A2  A3  A4  As  As  A^  ^CT  ^SO 

-3.78  -0.44  3.45  -1.63  1.68  -2.23  0.025  38  -13 


*  The  Ai  parameters  are  in  units  of  /l/2me  for  i=l. . .  6  and  e^/2  for  Aj,  Acr  and  Ago  are  in  meV. 

Spin  orbit  splitting 
Aso/3  A2  As 
-4.33  -4.53  -3.05 


the  g-factor  directly  from  our  first-principles  band  structures  using  perturbation  theory.  We  hope  to 
finish  this  work  during  the  requested  nocost  extension. 

In  summary,  the  controversy  of  the  valence  band  ordering  revolves  around  the  sign  of  the  spin- 
orbit  splitting,  which  depends  sensitively  on  the  Zn  d-band  position.  Our  calculations  including 
considerations  beyond  LDA  strongly  favors  a  negative  spin  orbit  splitting.  This  disagrees  with  the 
current  interpretation  in  the  literature  of  the  photoluminescence  Zeeman  splittings.  However,  the 
latter  is  based  essentially  on  an  assumption  about  the  sign  of  the  hole  g-factor  which  for  free  holes 
has  not  been  determined  independently.  An  alternative  explanation  reconciling  the  theory  with  the 
experiment  is  possible  if  the  hole  g-factor  of  this  particular  valence  band  is  negative.  Work  is  in 
progress  to  determine  the  sign  of  the  hole  g-factor.  Another  possible  reconciliation  between  theory 
and  experiment  is  that  the  assigment  as  a  Ps  and  Fe  exciton  would  be  correct  but  that  this  exciton 
ordering  cannot  be  directly  equated  to  the  valence  band  ordering  because  of  the  intermixing  of  excitons 
derived  from  different  valence  band  states.  Such  intermixing  is  to  be  expected  because  of  the  small 
band  splitting  compared  to  the  the  large  exciton  binding  energy.  This  aspect  of  the  problem  requires 
further  investigation.  An  unraveling  of  this  problem  in  conjuction  with  Reynold’s  data  will  further 
provide  information  on  the  exciton  exchange  coupling  paramter. 

Finally,  we  note  that  similar  questions  recently  have  appeared  for  GaN.  The  fine  splittins  due 
to  the  exchange  interaction  of  excitons  have  recently  been  observed  in  GaN.  In  GaN,  as  in  ZnO, 
the  binding  energy  of  the  excitons  is  comparable  and  exceeds  the  valence  band  splittings,  hence  an 
explicit  inclusion  of  the  exciton  intermixing  effects  must  be  included  before  conclusions  on  the  valence 
band  splittings  can  be  drawn.  This  may  in  part  help  to  resolve  the  discrepancies  between  current 
LDA  theory  for  the  crystal  field  splitting  and  the  values  deduced  from  experiment.  We  plan  also  to 
reinvestiage  the  effects  of  the  Ga  d-bands  in  view  of  what  we  have  learned  from  the  study  of  ZnO. 

As  part  of  this  work,  we  will  have  obtain  the  valence  band  effective  Hamiltonian  mass  parameters, 
in  ZnO.  Our  present  best  set  of  parameters  for  the  Rashba-Sheika-Pikus  Hamiltonian,  described  e.g. 
in  ref.  [18]  are  given  in  Table  2. 

3  Mga;Zni_;cO  alloys. 

Our  first  model  of  the  band  gap  bowing  in  the  Mga:Zni_j;0  alloy  system  is  based  on  calculations 
for  zincblende  based  structures.  In  particular,  we  used  a  50%  ZnO,  MgO  structure  and  the  end 
compounds.  In  this  model,  the  residual  strain  in  the  alloy  is  taken  into  account  by  considering  the 
expected  bondlengths  of  Zn-0  and  Mg-0  based  on  a  Keating  model  for  the  dilute  limits.  This  method 
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was  developed  in  earlier  work  under  this  grant  and  is  explained  in  detail  in  Ref.  [19].  Prom  a  statistical 
treatment,  the  band  gap  bowing  is  then  estimated  for  zincblende  ZnMgO.  It  is  found  to  be  0.56  eV. 
Determining  the  band  gaps  of  the  end  compounds  is  not  trivial  because  MgO  has  naturally  rocksalt 
structme.  We  calculated  the  band  structure  in  the  rocksalt,  wurtzite,  relaxed  wurtzite  (see  below) 
and  zincblende  structures.  We  calculated  ZnO  in  the  wurtzite  and  zincblende  structures.  All  of 
these  calculations  were  in  the  LDA.  Prom  the  known  experimental  gap  in  rocksalt  and  assuming  the 
correction  to  be  independent  of  structure,  we  then  obtain  an  estimated  band  gap  of  wurtzitic  IVIgO. 
Along  with  the  experimental  ZnO  gap  and  the  bowing  coefficient,  this  provides  the  calculated  data 
for  the  Mga;Zni_a;0  alloys  shown  in  Pig.  3.  In  spite  of  the  numerous  assumptions  this  calculation  is 
in  excellent  agreement  with  the  data  by  Ohtomo  et  al.  [2]. 

Although  these  alloys  have  been  successfully  grownt  up  to  about  40  %,  they  were  found  in  annealing 
studies  to  be  only  stable  towards  phase  separation  up  to  about  15  %.[20]  Also,  from  the  theory  point  of 
view,  many  precautions  are  in  order.  To  further  investigate  the  stability  of  the  alloys  in  wurtzitic  vs. 
rocksalt  structure  and  the  effects  of  local  relaxation  on  the  bowing,  we  started  with  an  investigation 
of  the  stability  of  the  end  compound  MgO. 

To  our  surprise,  what  we  found  is  that  the  MgO  wurtzite  structure  is  not  stable  at  all.  In  fact, 
restricting  the  structure  to  wurtzite  symmetry,  the  u  parameter  which  gives  the  ratio  of  the  bond 
length  in  the  c-direction  to  the  c  lattice  constant,  was  found  to  relax  to  a  value  of  1/2.  What  this 
implies  is  that  the  structure  obtains  a  additional  mirror  symmetry  and  degenerates  into  the  h-BN 
structure.  This  striu  'ura  i?.  stSI  hexago..?!  but  now  the  O  atoms  above  and  below  Mg  are  at  equal 
distance  and  the  bn  ..ling  of  the  Iexa.gc  :k:  network  is  wiped  out.  The  hexagonal  honeycomb  layers 
become  completely  flat.  Typically  in  h-BN  the  interplanar  distance  is  much  larger  than  the  in-plane 
distance,  resulting  in  effectively  3-fold  coordination  typicaxlly  of  7r-bonding  in  graphite.  However, 
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here  the  bonds  in  the  c-direction  are  about  the  same  as  in  the  plane.  This  means  that  Mg  is  roughly 
5-fold  coordinated  instead  of  4-fold.  This  can  be  understood  by  the  preference  of  ionic  materials  for 
high  coordination.  Within  the  restrictions  of  hexagonal  symmetry  this  is  the  best  MgO  can  do.  It 
was  then  discovered  that  a  further  in-plane  deformation  changes  the  structure  to  rocksalt.  Changing 
the  6/a  ratio  where  6  refers  to  the  distance  \/3o  in  the  [1200]  direction  to  a  will  turn  the  rectangular 
2x1  projected  cell  into  a  square  projected  cell  of  the  rocksalt.  The  energy  as  a  function  of  the 
transformations  is  shown  in  Fig.  4. 

The  structures  are  shown  in  Fig.  5. 

We  see  that  the  h-BN  structure  has  a  metastable  region  only.  This  discovery  of  a  continuous  path 
of  tranformation  between  wurtzite  and  rocksalt  prompted  us  to  investigate  the  same  transformation  in 
GaN.  In  fact,  in  GaN,  and  other  nitrides,  a  transformation  from  wurtzite  to  rocksalt  has  been  observed 
at  high  pressures.  The  energy  landscape  for  this  transformation  with  a  plausible  transformation  path 
is  shown  in  Fig.  6  as  a  contour  plot.  It  shows  that  the  rocksalt  structure  has  to  overcome  a  barrier  to 
transform  back  to  wurtzite.  Note  again  that  in  MgO  the  situation  is  inverted:  rocksalt  is  the  stable 
structure,  wurtzite  is  unstable. 

This  raises  serious  questions  on  how  the  structure  will  locally  relax  aroudn  Mg  atoms  in  ZnMgO 
alloys.  This  question  is  now  under  further  investigation. 

Further  details  on  the  rocksalt  to  wurtzite  transformation  will  be  given  in  a  paper  that  is  currently 
in  preparation. 

Another  important  parameter  for  the  use  of  Mga;Zni_3;0  alloys  as  barrier  material  along  with  ZnO 
quantum  wells  is  the  band  offset.  We  found  that  the  band  offset  is  almost  entirely  in  the  conduction 
band.  This  was  conifrmed  by  data  of  Ohtomo  et  al.  [20]  who  found  a  ratio  of  AEy/SEc  ~  1/9.  This 
means  that  electrons  will  have  good  confinement  but  holes  will  not.  In  view  of  this  unfavorable  result, 
we  have  also  investigated  some  other  poissible  barrier  materials  and  their  band  offsets  to  ZnO.  Amont 
these  an  interesting  candidate  is  LiGa02,  a  materials  closely  related  in  structure  to  ZnO.  All  these 
band  offset  calculations  were  perforemed  in  the  simple  dielectric  midgap  model  and  require  refinement 
to  take  in  to  consideration  the  specific  strains  in  lattice  mismatched  situation  and  the  specific  interface 
orientations.  However,  they  are  known  to  provide  at  least  initial  guidance  for  other  semiconductor 
materials.  Our  calculated  band  offsets  are  shown  in  Fig.  7 
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WZ  A-BN  RS 


Figure  4:  Total  energy  of  MgO  under  the  wurtzite  to  h-BN  and  h-BN  to  rocksalt  continous  transfor¬ 
mation.  The  volume  and  the  internal  structural  parameters  are  relaxed  for  each  c/a  or  b/a  ratio. 


Figure  5:  From  left  to  right,  ideal  wurtzite  structure,  h-BN  structure  (or  wurtzite  with  u  =  1/2)  with 
bond  length  ratios  corresponding  to  the  metastable  structure  of  MgO,  and  rocksalt  structure. 
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Figure  6:  Total  energy  of  GaN  as  function  of  c/a  and  b/a  parameters,  as  defined  in  text  with  relaxed 
volume  and  internal  parameters. 


AIN 


Figure  7\  Band-cusets  between  ZnO  and  related  materials.  The  dashed  line  gives  the  dielectric  midgap 
point  with  respect  to  which  the  band  structures  are  aligned. 
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4  Publications  and  reports. 

The  following  papers  were  published  or  are  in  preparation. 

1.  W.  R.  L.  Lambrecht,  S.  Limpijumnong,  and  B.  Segall,  Theoretical  studies  of  ZnO  and  related 
MgiZni-iO  alloys,  in  GaN  and  Related  Alloys,  edited  by  Stephen  J.  Pearton,  Chhping  Kuo, 
Alan  F.  Wright,  Takeshi  Uenoyama,  Mater.  Res.  Soc.  Symp.  Vol.  537,  p.  G6.8,  also  appearing 
in  MRS  Internet  J.  Nitride  Semicond.  Rese  4S1  G6.8  (1999). 

2.  H.  MorkoQ,  R.  Congolani,  W.  Lambrecht,  B.  Gil,  H-X.  Jiang,  J.  Lin,  D.  Pavlidis,  and  K.  Shenai, 
Material  Properties  of  GaN  in  the  Context  of  Electron  Devices,  in  GaN  and  Related  Alloys, 
edited  by  Stephen  J.  Pearton,  Chhping  Kuo,  Alan  F.  Wright,  Takeshi  Uenoyama,  Mater.  Res. 
Soc.  Symp.  Vol.  537,  p.  G1.21,  also  appearing  in  MRS  Internet  J.  Nitride  Semicond.  Rese  4S1 
Gl.2.1  (1999). 

3.  Theoretical  study  of  wurtzite  to  rocksalt  tranformation  in  MgO  and  GaN,  in  preparation  for 
Phys.  Rev.  B 

The  following  presentations  were  given  at  conferences. 

1.  MRS  Fall  meeting  1999,  poster  G6.8  in  Symposium  S4.1 

2.  APS  March  Meeting 

3.  First  International  ZnO  workshop  in  Dayton,  OH,  October  7-8,  1999 

Other  activities  during  this  period.  The  PI  during  this  period  also  was  co-author  on  a  paper  with 
Morkog  et  al.  relating  to  the  piezoelectric  and  spontaneous  polarization  effects  in  the  nitrides  and 
their  impact  on  devices. 

He  also  consulted  on  a  paper  on  the  electronic  transport  in  ZnO  through  a  collaboration  with  P. 
P.  Ruden  and  J.  D.  Albrecht  at  the  University  of  Minnesota.  The  important  result  of  this  work  is 
that  ZnO  has  a  higher  drift  velocities  at  high  fields  than  GaN  and  that  this  is  due  to  the  difference  in 
energy  from  the  conduction  band  to  the  first  secondary  value  being  larger  in  ZnO  than  in  GaN.  The 
work  was  mainly  carried  out  at  Univ.  Minnesota  but  was  based  on  the  band  structure  calculations  in 
our  group  at  CWRU. 

This  resulted  in  a  paper  submitted  to  J.  Appl.  Phys.  “High  field  electron  transport  properties  of 
bulk  ZnO,”  J.  D.  Albrecht,  P.  P.  Ruden,  S.  Limpijumnong,  W.  R.  L.  Lambrecht,  and  K.  F.  Brennan 

The  graduate  student  Sukit  Limpijumnong  involved  in  the  work  on  ZnO  presented  his  Ph.D  thesis 
in  June  1999  with  work  on  SiC  polytypes.  The  present  grant  has  contributed  to  his  support  during 
this  period. 
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